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Abstract. We present an analysis of COMPTEL obser- 
vations made between November 1991 and May 1994 
of 2CG 135+01, a bright 7-ray source located near the 
Galactic plane. At energies above 1 MeV, an excess con- 
sistent with the position of 2CG 135+01 is detected in the 
sum of the observations, at flux levels which are a factor 
of 10 — 100 below those published in the past. The de- 
tection significance of this excess, when the possible pres- 
ence of underlying Galactic diffuse emission is neglected, 
is 6.6(7 for 3 degrees of freedom. The differential photon 
spectrum in the 1-30 MeV energy range can be described 
by a power law with a spectral index of 1.95^q'3. Due 
to the uncertainties involved in modelling the Galactic- 
disk diffuse emission underneath the source, the absolute 
flux levels must be considered uncertain by a factor of 
two. They are consistent with the extrapolation of the 
time-averaged spectrum of 2CG 135+01 measured with 
EGRET, thereby strengthening the identification. No sig- 
nificant temporal correlation between the 7-ray emission 
and the monitored radio emission of the possible coun- 
terpart radio source GT 0236+610 (showing a 26.5 day 
modulation) is found. 

Key words: Gamma rays: observations - Stars: individ- 
ual: 2CG 135+01 - Stars: individual: GT 0236+610 
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1. Introduction 

With a flux above 100 MeV of 1.0 x 10"*^ photons cm'^ 
s~^, the 7-ray source 2CG 135+01 is one of the bright- 
est unidentified high-energy Galactic-plane sources in the 
second COS-B catalogue (Hermsen et al. 1977; Swanen- 
burg et al. 1981). 2CG 135+01 attracted attention when 
its position was found to be consistent with that of the 
radio source GT 0236+610, which exhibits strong radio 
outbursts of a non-thermal character (Gregory & Tay- 
lor 1978). Long-term radio observations of GT 0236+610 
subsequently revealed a 26.496-day periodicity in these 
outbursts (Taylor & Gregory 1982, 1984), of which the 
amplitude is possibly modulated by a 4-year period (Gre- 
gory et al. 1989; Paredes et al. 1990). The decay of the 
radio fiux during outbursts is reminiscent of synchrotron 
emission from an expanding cloud of relativistic electrons 
(Taylor & Gregory 1984). In a high-resolution map ob- 
tained two days post-outburst from VLBI measurements 
at 6 cm, the object appears to consist of two compo- 
nents separated by 3.1 x 10"'^^(I?/2.3 kpc) cm (Massi et 
al. 1993). Optically, GT 0236+610 has been identified 
with LS I +61°303, which has a spectrum typical of an 
early-type B star (Hutchings & Crampton 1981; Maraschi 
et al. 1981; Paredes & Figueras 1986) and exhibits the 
near-infrared excess commonly observed for Be-type stars 
(DAmico et al. 1987; Hunt et al. 1994). Recent modeUing 
of the JHK-band light curves of LS I +61°303 has shown 
that the onsets of the radio outbursts roughly coincide 
with the inferred periastron passage (Martf & Paredes 
1995). The source has also been detected in X-rays (Big- 
nami et al. 1981; Goldoni & Mereghetti 1995). The X-ray 
spectrum is rather hard compared to those of normal B 
stars, but is similar to spectra observed for Be/X-ray bi- 
naries. 
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In recent observations made with EGRET on board 
the Compton Gamma- Ray Observatory (CGRO), a som-ce 
designated 2EG J0241-I-6119 was detected on several oc- 
casions (Thompson et al. 1995). This source, whose po- 
sition is consistent with that of 2CG 135-1-01 and which 
is roughly as bright, is generally assumed to be the COS- 
B source. The improved position of the 7-ray source es- 
tabhshed by the EGRET detection of 2EG J0241+6119 
(a 95% error radius of 33') is consistent with that of 
GT 0236-^610 (located at {l,b) = (135.68°, 1.09°)) and 
has strengthened the proposed identification with the ra- 
dio source. In addition, the nearby QSO 4U 024H-61 (lo- 
cated at {l,b) = (135.64°, 2.43°)), which was also con- 
tained in the relatively large COS-B error region, is now 
firmly rejected (von Montigny et al. 1993). 

A detection of the 26.496 day periodicity in 7-rays 
would remove any remaining doubt about the identifica- 
tion with GT 0236-f 610, but neither COS-B nor EGRET 
has yet detected such flux variations. Note that the 
26.496-day periodicity inferred from radio observations of 
GT 0236-1-610 has possibly been detected at infrared and 
optical wavelengths as well (Mendelson & Mazeh 1989; 
Paredes et al. 1994), with smaller flux variations towards 
shorter wavelengths. At UV wavelengths and in X-rays, 
where the positional accuracy is good enough to allow a 
confident identification with GT 0236-1-610, orbital-phase 
related flux variations have not been detected up to re- 
cently (Bignami et al. 1981; Howarth 1983; Goldoni & 
Mereghetti 1995). The factor of ~ 10 increase of the X-ray 
flux in the light curve measured with ROSAT during the 
monitoring of one orbital cycle (Taylor et al. 1996) may 
be the first evidence of orbital- modulated X-ray emission. 

Despite the firm detections of GT 0236+610 up to 
X-rays and of 2CG 135+01 above 100 MeV, the situa- 
tion around 1 MeV is less clear. When we compare the 
spectral index measured for GT 0236+610 with ROSAT 
{a w 1.1: Goldoni & Mereghetti 1995) with that mea- 
sured for 2CG 135+01/2EG J0241+6119 with EGRET 
(a « 2.21: Fierro 1995), it is evident that a spectral break 
must occur in between. This naturally also holds if the two 
sources are not related. A simple extrapolation of the mea- 
sured power laws to the MeV energy range predicts fluxes 
in the range of 10^^ — 10^^ photons cm~^ s~^ MeV~^, 
which is close to the detection threshold of COMPTEL. 
A preliminary analysis of COMPTEL observations of the 
2CG 135+01 field, based on three observations made be- 
tween November 1991 and August 1992, did not reveal ev- 
idence for such a strong source but showed the presence of 
a weak excess consistent with the position of 2CG 135+01 
(van Dijk et al. 1994). 

2. Instrument and Observations 

COMPTEL (Schonfelder et al. 1993) is a 7-ray telescope 
on board CGRO, sensitive in the 0.75-30 MeV range. It 



has a wide field of view of 1 steradian, a source location 
accuracy of 1° and an energy resolution ranging from 
10% at 1 MeV to ^ 5% at 6 MeV. Incoming photons 
first Compton scatter on an electron in the upper layer 
of detectors (Dl) through an angle (pg after which they 
are completely or partially absorbed in the lower detector 
layer (D2). 

The COMPTEL data consist of events which are de- 
fined as coincident interactions in the two detector layers. 
These events are binned in a 3-dimensional data space, 
of which the first two dimensions denote the coordinates 
X and of the direction of the scattered photon. The 
third dimension, the calculated scatter angle if, follows 
from the well-known Compton formula and is equal to the 
true scatter angle ipg only if the electron and photon are 
completely absorbed in Dl and D2 respectively. 

Given a celestial intensity distribution, the expected 
distribution of events in the 3D data space is the convo- 
lution of the intensity distribution with the instrument 
response function, which is the product of a 3D point- 
spread function (PSF) and a total exposure. The latter 
has been subdivided into a 3D geometry function and a 
2D exposure function for practical reasons. Throughout 
this paper we used PSFs obtained from Monte Carlo sim- 
ulations, taking into account the full spherical response 
representation. 

Table |l] gives a log of the observations which were 
used for the analysis presented here. Unfortunately, during 
most of these observations, the angular distance between 
the source and the pointing direction was not favorable. 



Obs 


Time 


Phase 


zn 


15 


1991, Nov. 28-Dec. 12 


0.096-0.631 


22.2 


31 


1992, Jun. ll"Jun. 25 


0.499-0.023 


29.3 


34 


1992, Jul. 16-Aug. 06 


0.820-0.610 


27.2 


211 


1993, Feb. 25-Mar. 09 


0.272-0.725 


11.4 


319.0 


1994, Mar. 01-Mar. 08 


0.197-0.464 


28.1 


319.5 


1994, Mar. 15-Mar. 22 


0.727-0.990 


27.0 


325 


1994, Apr. 26-May 10 


0.311-0.840 


15.2 



Table 1. List of CGRO observation periods during which 
2CG 135+01 was within 30° of the pointing direction. The 4 
columns give: 1) the observation number in CGRO notation; 2) 
the start and end date of the observation; 3) the corresponding 
radio phase interval; 4) the angular distance of 2CG 135+01 
to the pointing direction. 



3. Analysis 

The COMPTEL data were analysed in the standard en- 
ergy ranges 0.75-1, 1-3, 3-10 and 10-30 MeV, as well as 
in a combined 1-30 MeV energy range. Significance and 
flux maps were produced using the maximum-likelihood 
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Fig. 1. This figure shows the (x, '!/')"distribution of events in 
a data space containing simulated point sources at the posi- 
tion of the Crab and at the refined location of 2CG 135-1-01 / 
2EG J0241-(-6119, subject to a selection on the derived scat- 
ter angle of [26°, 28°]. In this representation, point sources 
show up as ring-shaped distributions centered on the source 
positions with a radius of ip°. 

ratio (MLR) method. For a description of the application 
of this method to COMPTEL data, the reader is referred 
to de Boer et aL (1992). The Icr errors quoted in this paper 
are statistical only, i.e. they do not include the estimated 
systematic 30% flux uncertainty. 

3.1. Model fitting and the instrumental background model 

To account for flux contributions from nearby point 
sources and the underlying Galactic diffuse emission, the 
MLR analysis allows for the inclusion of data-space models 
in addition to an instrumental background model. Such 
data-space models are created by folding the celestial 
intensity distribution, a simple delta function for point 
sources, through the instrument response. 

The basic principle of the instrumental background 
model used here is described by Bloemen et al. (1994). 
It is created by fitting (^-templates, obtained from the 
scaled geometry function, to the event data space itself. 
This background model has recently been improved upon 
by iteratively correcting the (^-templates for contributions 
from fitted models of detected celestial sources. The modi- 
fication becomes important when there are strong sources 
(e.g., the Crab or large-scale diffuse emission) in the field 
of view. In the first iteration of the model fitting, no cor- 
rection is performed. If a positive signal is found for one or 
more models, a second iteration can be performed during 
which the (^-templates are corrected for the (^-dependent 
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Fig. 2. Location contours (la, 2a, 3a) for the 1-30 MeV excess 
for the sum of the observations listed in Table |l[ The positions 
of GT 0236-^610 and the QSO 4U 0241-^61 are denoted with 
and '-f ' respectively. 

model contributions. This procedure is repeated until the 
scale factors of the fitted models converge. 

The effective exposure function for the observations 
listed in Table |l| covers a large part of the sky and extends 
beyond I = 200°. We therefore must take into account the 
non-negligible influence of the events originating from the 
Crab 2.5 x 10'' in the 1-30 MeV range compared to ~ 
7 X 10^ for 2CG 135-1-01) on the background model by in- 
cluding a data-space model for the Crab. For this set of ob- 
servations, the situation is particularly confusing because 
the signals of the Crab and of a source at the position of 
2CG 135+01/2EG J024H-6119 are at a maximum in the 
same region of the data space (the region where the geom- 
etry function peaks, see Fig. |l|). When models for both the 
Crab and for 2CG 135-f 01/2EG J0241-h6119 are included 
simultaneously, this sometimes causes the Crab flux to be 
underestimated and the 2CG 135-h01/2EG J0241-h6119 
flux to be overestimated in the flrst iteration of the back- 
ground modelling. In general, 4 iterations arc needed to 
correct for this effect. 

4. Results 

4-.1. Spatial analysis 

When all the observations in Table |^ are summed, an ex- 
cess consistent with the position of 2CG 135-1-01 is found 
at energies above 1 MeV. Below 1 MeV an upper limit 
is obtained. Fig. || shows the location contours for the 1- 
30 MeV energy range; the maximum of the excess has a 
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MLR value of 50.9. This corresponds to a point-source 
detection significance of 6.6ct for 3 d.o.f. (neglecting the 
possible presence of underlying Galactic diffuse emission), 
while at the position of 2CG 135-hOl it is 6.7(t (1 d.o.f.). 
In the standard energy ranges 1-3 MeV, 3-10 MeV and 
10-30 MeV, the peak MLR values are 24.4, 18.8 and 10.3 
respectively. In the individual observation periods the sig- 
nal is either weak or not detected at all. 

In Fig. the significance maps for the standard en- 
ergy ranges above 1 MeV are shown. In each of these maps, 
an excess consistent with the position of 2CG 135-t-Ol 
can be observed. The significance in the highest energy 
range is marginal (2.6(t), while below 3 MeV the ex- 
cess is broader than expected for a single point source. 
When a data-space model for a source at the position of 
2CG 135-f 01/2EG J0241-f 6119 is included, we obtain the 
residual maps as shown in Fig. ^jb. None of the remain- 
ing excesses by itself is formally significant at the > 3cr 
level for a free search, but the number of 2a — 3cr ex- 
cesses yield it unlikely that all are just statistical fluctu- 
ations. One of the possible celestial components, tha t of 
the Galactic diffuse emission, is discussed in Section 4.4. 



we show the results for the extreme values of index a as 
well as for a = 2. Also listed are the expected ratios for 
power-law spectra with these indices. 



A possible contribution from point sources, however, can- 
not be excluded. Most noticeable is the weak emission 
around (/, b) = (134.6°, -^7.0°) = Pi in the 1-3 MeV map 
in Fig. y, which is completely due to a strong (3.9cr) excess 
observed during Obs. 325. It is not seen in any other en- 
ergy range nor in any other COMPTEL observation (see 
also Bhattacharya & Owens 1994). After a correction for 
a minimum number of trials (4 energy intervals, 7 obser- 
vations), its significance in Obs. 325 is only slightly larger 
than 3ct, which does not confidently exclude a statistical 
background fluctuation as origin. However, located only 6° 
away from 2CG 135+01/2EG J0241+6119 (but not con- 
sistent with the position of the latter at the > 3(T level), 
the apparent source Pi has a non-negligible influence on 
the flux derived for the COMPTEL source in the mosaic 
of the observations and thus has to be accounted for. 

4-2. Spectral analysis 

The analysis in the previous section was performed with 
PSFs for an assumed power-law photon spectrum 

with a = 2. Since the fluxes and likelihoods depend only 
slightly on the assumed input spectrum, this was not a bad 
approximation. In order to find the power-law index that 
fits the data best, the MLR analysis has been repeated 
using PSFs with various values for a until the ratios of the 
fiuxes measured in different energy ranges are consistent 
with the assumed input spectrum. 

Such flux ratios were calculated using PSFs with in- 
dices ranging from 1.4 to 2.8. In the model fltting we in- 
cluded a data-space model for the Crab in each energy 
range, while for the 1-3 MeV range we also included a 
model for the nearby excess PI to correct for its influence 
on the flux derived for the COMPTEL source. In Table | 



Index 


Ra (1-3/3-10) 


Rb (3-10/10-30) 


Rc (1-3/10-30) 


1.4 


2.3 ±0.8 


4.1 ± 1.6 


9.4 ±3.8 




1.44 


1.74 


2.51 


2.0 


2.2 ±0.8 


4.1 ± 1.6 


9.1 ±3.8 




2.86 


3.50 


10.00 


2.8 


2.1 ±0.8 


4.2 ± 1.6 


8.7 ±3.8 




7.03 


8.98 


63.10 



Table 2. Measured (first line of each entry) and expected (sec- 
ond line of each entry) hardness ratios for power-law input 
spectra with different spectral indices. These hardness ratios 
are defined as the ratio of the fluxes in the specified energy 
ranges (units [photons cm~^ s~^]). 



E [MeV] 


Flux [photons cm ^ s ^] 




no diff. models 


incl. diff. models 


0.75-1 


< 4.8 X 10"'' 


< 4.8 X 10"'' 


1-3 


(11.2 ± 3.2) X 10"^ 


(6.0 ±3.0) X 10"^ 


3-10 


(5.0 ± 1.2) X 10"^ 


(2.2 ± 1.2) X 10"^ 


10-30 


(1.2 ±0.4) X 10"^ 


(0.3 ±0.4) X 10-5 



Table 3. Time-averaged fluxes with la errors and the 2a upper 
limits for the COMPTEL source consistent with the position 
of 2CG 135±01, for the sum of observations given in Table |l| 
These results were derived including a model for the Crab in 
all energy ranges. In the 1-3 MeV range a point-source model 
at {l,b) = (134.6° , 7.0°) was included as well (see text). The 
second column shows the results obtained disregarding the pos- 
sible presence of underlying diffuse emission, the third column 
shows the results when models for the Galactic diffuse emission 



are incorporated (see Section 4.4) 



Although the measured ratios are not very sensitive 
to the assumed input spectrum, the expected ratios dif- 
fer considerably. The indices that flt the COMPTEL data 
best are ao — I.SIq;!^ and ao — 2.15j;o;4 for the Ra and 
Rb hardness ratios respectively. They are consistent with 
a single power law F = I{E/1 MeV)""" across the entire 
COMPTEL energy range with ao = 1.95to 3 ^^'^ ampli- 
tude / = 1.75x10"'' photons cm"^ s"^ MeV"^ Adopting 
this spectral index, the time-averaged COMPTEL fluxes 
and upper limits are given in the second column of Table |[ 
The upper limit in the 0.75-1 MeV range is inconsistent 
with the power law at the 2.4ct level, which suggests that 
there is a spectral break in the COMPTEL energy do- 
main. Note that at 1 MeV, the COMPTEL flux levels are 
a factor of 10 — 100 below those reported in the past for 
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Fig. 3. MLR maps for the region around 2CG 135+01, centered on {l,b) — (135.5° , 1.5°) with a radius of 20°, longitude I 
increasing towards the left. The horizontal and vertical lines are drawn at fo = 0° and I — 135° respectively and the position 
of 2CG 135+01/2EG J0241+6119 is denoted by the 'x' symbols. Contour levels are spaced as 4, 9, 16 and 25. a) including a 
model for the Crab; b) including models for both the Crab and for 2CG 135+01/2EG J0241+6119. 



the emission from this celestial region (Coe et al. 1978; 
Perotti et al. 1980). 



4-3. Temporal analysis 

We performed a search for flux variations correlated with 
the orbital period of GT 0236+610. For this the data 
for each energy range were binned and analysed in 5 
phase bins of width 0.2, assuming an orbital period of 
26.496±0.008 days with phase at JD 2,443,366.775 (Tay- 
lor & Gregory 1984). In the phase bins that contained a 
contribution from Obs. 325, we again included a point- 
source model for the Pi excess observed in that observa- 
tion in the 1-3 MeV range (see Section 4.1). The radio 



lightcurves differ significantly from cycle to cycle and the 
phase of the peak outburst is seen to vary by more than 
^ 0.5. Therefore, this approach may fail to detect correla- 
tions of the weak 7-ray flux with the orbital period if the 
7-ray flux variations occur, like the radio-flux variations, 
at irregular phases. Note that the error on the orbital pe- 
riod corresponds to a total relative phase shift during the 
observations given in Table of only 1%, but that the ab- 
solute phase may be wrong by ~ 6%, corresponding to 
^ 1.5 days. 



The results of this analysis are shown in Fig. ^. In each 
energy range, a x^-fit with a constant distribution has 
been used to determine if there is evidence for orbital- 
phase related flux variations. The values thus obtained 
are 5.49, 6.10 and 3.38 for 4 d.o.f., indicating that the 
COMPTEL data do not show evidence for flux variations 
correlated with the orbital phase. Due to the large statis- 
tical flux errors, however, and the above-mentioned pos- 
sibility of peak fluxes occurring at irregular phases, 7-ray 
variability similar to that observed in X-rays with ROSAT 
(Taylor et al. 1996) cannot be excluded. 

4-.4- The Galactic diffuse emission 

So far in the analysis we have not taken into account the 
possible contribution of diffuse emission from the Galac- 
tic disk. If this component, which is dominant at energies 
above 100 MeV, is indeed not neghgible at COMPTEL 
energies in this region of the sky, then the flux of the 
COMPTEL source could be significantly lower than de- 
rived above. Furthermore, Monte Carlo simulations have 
shown that weak diffuse emission from a Galactic-disk dis- 
tribution may result in several weak point-source-like ex- 
cesses at low latitudes, similar to those observed in Fig. |[ 
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Fig. 4. Flux [ photons cm~^ MeV~^] as a function of or- 
bital phase for three energy ranges. The dashed hue shows the 
result from a x^"fit with a constant distribution. 



4.4.1. The diffuse models 

The classical interpretation of the Galactic diffuse emis- 
sion above 100 MeV is that it is dominated by 7r°-decay, re- 
sulting from nuclear interactions between cosmic-ray (CR) 
particles and nuclei of the interstellar gas (for an overview 
see Bloemen 1989). Below 100 MeV, the main contribu- 
tion is thought to come from bremsstrahlung losses of 
low-energy CR electrons interacting with the interstellar 
gas. In both cases, the celestial 7-ray distribution is ex- 
pected to correlate approximately with that of the inter- 
stellar gas, as mapped by HI and CO surveys (the latter 
component assumed to be a tracer of H2). Studies per- 
formed with SAS-2 and COS-B data (^ 100 MeV) have 
shown that such a correlation indeed exists (e.g., Fich- 
tel & Kniffen 1984; Bloemen 1989). It should be noted 
that the assumptions about the CR density distribution, 
e.g., whether or not it is coupled to matter, differ among 
the various studies. Another component (= IC) that is 
thought to contribute to the Galactic diffuse emission both 
below and above 100 MeV is inverse Compton radiation 
from high-energy CR electrons scattering on the interstel- 
lar radiation field. The relative contribution of this com- 
ponent near the Galactic plane probably does not exceed 
5%-10% at energies near 100 MeV (see Bloemen, 1989, 
for an overview). 

Recent analyses revealed that the COMPTEL data 
show evidence for a broader latitude distribution than 
those of the HI and CO models. It remains uncertain, 
however, whether this should indeed be attributed to IC 
emission. Furthermore, an emission component that cor- 
relates with the CO distribution appears to be lacking in 



the COMPTEL energy regime (Strong et al. 1996; Bloe- 
men et al. 1996), except perhaps for the 10-30 MeV energy 
range (van Dijk 1996). Also, it should be noted that part 
of the inferred diffuse intensities may be due to unresolved, 
weak point sources. Clearly, the problem of the Galactic 
diffuse emission at COMPTEL energies requires further 
study. Here we will use the current understanding to dis- 
cuss the pos sible imp lications for the quantities derived in 
Sections O and O. 



4.4.2. Implications for derived quantities 

In order to assess the possible influence of the Galactic dif- 
fuse emission underneath the COMPTEL source on the 
derived source flux, we repeated the MLR analysis pre- 
sented above including data-space models for the inter- 
stellar gas, using the same HI and CO surveys that were 
used in the COS-B studies (Bloemen 1989) and adapting 
the model for the IC distribution from Strong & Youssefi 
(1995). Firstly, we fitted only the HI and IC models, fixing 
the HI emissivities and the scaling factors for the IC mod- 
els at the values found by Strong et al. (1996). As a result, 
the fluxes in the 1-3 MeV and 3-10 MeV ranges are re- 
duced by a factor of 2 (third column in Table H) . Above 
10 MeV, the flux drops below la signiflcance. The photon 
spectral index derived from these reduced 1-3 MeV and 
3-10 MeV fluxes is consistent with that derived in Sec- 



tion 4.2, although not very constraining. Since the errors 



on the fluxes have remained roughly constant, the tem- 



poral analysis from Section 4.3 is expected to yield even 
lower x^-values and was therefore not repeated. 

Secondly, we investigated a worst-case scenario in 
which the CO model is added to the HI and IC mod- 
els with a fixed value for the CO-to-H2 conversion factor 
= 2.0. Although correlations of the 7-ray emission 
with the CO distribution have not been found for COMP- 
TEL data, fixing Xj at the value found above 100 MeV 
gives an estimate of the error made in case the CO com- 
ponent is present, but for some reason is not recognized. 
The inclusion of the CO model is found to result in a small 
decrease of the 1-3 MeV and 3-10 MeV fluxes inferred for 
the COMPTEL source of less than ~ 15%. 



The effect of including the models for the Galactic dif- 
fuse emission on the parameters derived for the COMP- 
TEL source can only be considered to be crudely esti- 
mated. It is therefore likely that the fluxes given in the 
second column of Table ^ are overestimated, but due to 
the uncertainties in modelling the Galactic diffuse emis- 
sion, the fluxes given in the third column of Table ^should 
be interpreted with caution as well. 



5. Discussion 
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Fig. 5. The phase-averaged fluxes (xE^) for the COMPTEL excess 'as is' (points labeled 'COMPTEL'; second column in 
Table and those corrected for the estimated contribution from the Galactic diffuse emission (points labeled 'COMPTEL(D)'; 
third column in Table compared with unabsorbed X-ray and 7-ray data for GT 0236-1-610 (lines labeled 'ROSAT'; Goldoni 
& Mereghetti 1995), 2CG 135-h01/2EG J0241-F6119 (lines labeled 'EGRET'; Fierro 1995) and 4U 0241-^61 (dashed lines). AU 
data shown, except for the COMPTEL data, are extrapolated power laws, from the energy range in which the measurements 
were made (thick lines), into the MeV range (thin lines). In each case, two power laws are shown, reflecting the uncertainties in 
the photon spectral indices. The X-ray data for 4U 0241+61 (dashed lines) were measured with Einstein (Turner et al. 1991), 
OSO 8 (WorraU et al. 1980) and EXOSAT (Turner & Pounds 1989). 



5.1. Identification of the COMPTEL source 

The celestial distributions of neither the HI nor the IC 
models can mimic a point-source-like feature at the po- 
sition of 2CG 135+01. So, despite the uncertainties in 
the absolute flux levels of the COMPTEL source due to 
the unknown contribution of the Galactic diffuse emis- 
sion, the mere presence of a source consistent with the 
position of 2CG 135+01 is evident. If we extrapolate the 
time-averaged high-energy 7-ray spectrum measured for 
2CG 135+01/2EG J0241+6119 (Fierro 1995; lines labeled 
'EGRET' in Fig. |), we find that the COMPTEL fluxes 
prior to the correction for the Galactic diffuse emission 
(points labeled 'COMPTEL' in Fig. |) are consistent with 
the extrapolations. When the models for the Galactic dif- 
fuse emission are included, the fluxes are still consistent, 
but appear to be somewhat low (points labeled 'COMP- 
TEL(D)'). Although this may indicate that the contri- 



bution of the diffuse emission is overestimated, a spectral 
break is expected to occur somewhere in this energy range. 

The extrapolations of the X-ray spectra measured for 
GT 0236+610 (lines labeled 'ROSAT' in Fig. |) and for 
the quasar 4U 0241+61 (the dashed lines) do not re- 
ject any of the two being the possible counterpart of the 
COMPTEL source. There are several arguments, however, 
which favour GT 0236+610 being the counterpart. Firstly, 
4U 0241+61 has only been detected up to 20 keV 
(Turner & Pounds, 1989) and is not seen with EGRET (> 
30 MeV). On the other hand, both ROSAT and EGRET 
have reported detections coincident with GT 0236+610, 
with location uncertanties that leave little room for other 
counterparts (Goldoni & Mereghetti 1995; von Montigny 
et al. 1993; Thompson et al. 1995). Secondly, the only ac- 
tive galactic nuclei (AGNs) that have been detected with 
COMPTEL (and EGRET) up to now are those which are 
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radio loud and have flat radio spectra. This is consistent 
with the findings from an analysis of EXOSAT ME data 
by Lawson et al. (1992), who found a steeper average 
2-10 keV spectrum for radio-quiet QSOs (photon index 
a = 1.9±0.11) than for radio-loud QSOs (a = 1.66±0.07). 
Although 4U 0241-1-61 has been shown to contain a com- 
pact radio core with an inverted radio spectrum, the ra- 
dio flux of ^ 0.20 Jy at 5 GHz is rather low (Apparao 
et al. 1978) and has only been detected at radio wave- 
lengths because it is relatively nearby (z = 0.0438). Un- 
less 4U 0241-1-61 behaves very differently from other radio- 
quiet QSOs at MeV energies, it is not likely that its spec- 
trum extends into the COMPTEL energy range. 

Apart from GT 0236-f 610 and 4U 0241-1-61, there are 
no other sources in the SIMBAD data base that can be re- 
garded potential counterparts, in the sense that there are 
no strong non-thermal X-ray or high-energy 7-ray emitters 
nearby. A possible contribution from the large molecular 
complex W3 at {l,h) « (134°, 1°), as part of the Galac- 
tic diffuse emission, was estimated to have a small impact 
(Section |4Z|). 



5.2. Comparison with theory 

The previous sections have shown that the identification 
of the COMPTEL source with 2CG 135+01 is hkely, but 
that the absolute flux levels are uncertain by a factor of 
two. Here we will address the possible implications for the 
main theoretical models for 2CG 135+01/GT 0236+610, 
assuming for now that the radio source and the 7-ray 
source are indeed one and the same. 

Based on the presence of a ~ 1.4 Mq compact object 
(Hutchings & Crampton 1981), two basic models exist to 
explain the variable radio emission and the possible high- 
energy emission of the Be/X-ray binary GT 0236+610. 
In the supercritical accretion (SA) model (Taylor & Gre- 
gory 1982, 1984; Taylor et al. 1992), a neutron star is 
in an eccentric orbit around the Be star, accreting from 
the slow high-density equatorial wind of the latter. For 
suitable chosen values of the wind parameters and the 
eccentricity, the radio outbursts arise from the energetic 
particles produced during the supercritical accretion that 
occurs near periastron passage. The four-year modulation 
and the varying phases of the peak radio fluxes may be 
explained by a varying speed of the equatorial wind of the 
Be star (Martf & Paredes 1995). However, it is not clear 
why the X-ray luminosity during outbursts is several or- 
ders of magnitudes below that for other Be/X-ray binaries 
(Taylor et al. 1996). 

Maraschi & Treves (1981), on the other hand, proposed 
the so called young-pulsar model, in which the outbursts 
and high-energy emission originate in the shock front that 
results from the collision of the stellar-disk wind and the 
relativistic pulsar wind. Such a model has also been used 
for the highly non-thermal 1-200 keV emission detected 



from the Be star/pulsar system PSR 1259-63 containing 
a 47 ms radio pulsar (Tavani, Arons & Kaspi 1994; Ta- 
vani et al. 1996). In the case of the PSR 1259-63 system, 
accretion onto the surface of the neutron star is inhib- 
ited by the strong radiation pressure from the pulsar's 
electromagnetic and relativistic particle wind. The high- 
energy emission is produced by synchrotron and inverse- 
Compton emission of shock accelerated particles. Shock- 
powered high-energy emission is of general importance in 
pulsar binaries and the 2CG 135+01 system may belong 
to the category of hidden pulsars (Tavani 1995). 

Both types of models involve a neutron star in an ec- 
centric orbit around the Be star and predict, to some de- 
gree, time variability of the 7-ray emission correlated with 
the orbital period. The COMPTEL data on the other 
hand are consistent with a constant flux, but the large 
statistical errors on the COMPTEL fluxes (Table |) do 
not allow for tight constraints on the variability. At 7-ray 
energies above 100 MeV, the chance that the flux from 
2CG 135+01/2EG J0241+6119 is constant between ob- 
servations is only 0.1% (Fierro 1995), while an increase 
by a factor of ~ 2 was observed within a single observa- 
tion (Kniffen et al., in preparation). In the young-pulsar 
model of emission, orbital modulation of the high energy 
emission depends on the poorly understood pitch angle 
evolution of radiating shock-accelerated particles of the 
pulsar wind. Upstream of the shock, the particles flow ra- 
dially away from the pulsar and therefore initially have 
a nonisotropic velocity distribution. For large values of 
the ratio of the electromagnetic and plasma flow energies 
upstream, the shocked particles acquire preferential direc- 
tions of motion (Gallant et al. 1992) and a dependence of 
the 7-ray emission on orbital phase might be observed. If, 
on the other hand, this ratio is low, the onset of the Weibel 
instability may isotropize the downstream flow and there- 
fore the 7-ray emission (Weibel 1959; Hoshino et al. 1992). 
Note that in both of these cases, the shock-powered emis- 
sion may be more dependent on the varying conditions of 
the pulsar cavity being compressed by the gaseous outflow 
from the companion star (Tavani, Arons & Kaspi 1994). 

The COMPTEL measurements presented here imply 
a luminosity in the 1-30 MeV energy range of 5.0 x 
1035(D/2 kpc)2 ergs s'^ (2.0 x lO^^ ergs s'^ if the mod- 
els for the diffuse emission are included). Together with 
the ROSAT and EGRET measurements, this implies a 
total luminosity above 0.1 keV of L « 1.5 x lO'^^ ergs s~^. 
For comparison, the inferred total energy in the rel- 
ativistic electrons responsible for the radio outbursts, 
which typically last several days, is only 2 x 10*^^ ergs 
(Taylor & Gregory 1984). A high-energy luminosity near 
\Qi(> gj-gg g-i jjjay be powered by shock emission of a rela- 
tivistic pulsar wind, as observed in the case of the plerionic 
Crab Nebula and for the binary Be star/pulsar system 
PSR B1259-63 (Tavani et al. 1996). Here, the deduced 
spindown luminosity of the underlying pulsar would be 
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of the order of 10^'^ ergs s~^ for a conversion efficiency 
into the 7-ray band of about 10% (similar to the con- 
version efficiencies observed for the Crab Nebula and for 
PSR B1259-63). The possible spectral break observed for 
2CG 135-t-Ol in the MeV energy range may, when more 
accurately determined, be used to constrain parameters 
such as the kinetic energy of the particles and the electro- 
magnetic energy flux upstream. 

6. Conclusions 

We have detected a COMPTEL source consistent with the 
position of 2CG 135+01/2EG J0241+6119 in the energy 
range 1-30 McV. The flux at 1 McV is roughly two orders 
of magnitude below the values reported in the past for a 
source in this celestial region. The COMPTEL spectral 
shape is consistent with the spectral break expected from 
extrapolations of the ROSAT and EGRET observations 
of GT 0236+610 and 2CG 135+01/2EG J0241+6119 re- 
spectively. 7-Ray observations with higher spatial resolu- 
tion, or showing correlated variability, are necessary to 
prove unambigously the association between the radio 
source GT 0236+610 and the 7-ray source 2CG 135+01. 
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